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Abstract. In this paper, the problem of the identification of the symmetry class of a given tensor 
is asked. Contrary to classical approaches which are based on the spectral properties of the linear 
operator describing the elasticity, our setting is based on the invariants of the irreducible tensors 
appearing in the harmonic decomposition of the elasticity tensor [19] . To that aim we first introduce 
. a geometrical description of Ela, the space of elasticity tensors. This framework is used to derive 

' invariant-based conditions that characterize symmetry classes. For low order symmetry classes, such 

conditions are given on a triplet of quadratic forms extracted from the harmonic decomposition 
of C, meanwhile for higher-order classes conditions are provided in terms of elements of H*, the 
h-^ higher irreducible space in the decomposition of Ela. Proceeding in such a way some well known 

, conditions appearing in the Mehrabadi-Cowin theorem for the existence of a symmetry plane [16] are 

' retrieved, and a set of algebraic relations on H* characterizing the orthotropic ([D2]), trigonal ([D3]), 

tetragonal ([D4]), transverse isotropic ([S0(2)]) and cubic ([O]) symmetry classes are provided. Using 
' ' , a genericity assumption on the elasticity tensor under study, an algorithm to identify the symmetry 

' class of a large set of tensors is finally provided. 
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Introduction 



Physical motivation. In the theory of hnear elasticity, the stress tensor a and the strain tensor 
£ are related, at a fixed temperature, by the Hooke 's law 



If the medium under study is homogeneous, its elastic behavior is fully characterized by a 4th-order 
elasticity tensor C. The infinitesimal strain tensor e is defined as the symmetric displacement 
gradient 



and, assuming that the material is not subjected to volumic couple, the associated stress tensor is 
the classical symmetric Cauchy's one. Hence, both a and e belong to the 6-D vector space S'^(M^), 
where 5^ denotes the symmetric tensor product. As a consequence, the elasticity tensor is endowed 
with minor symmetries: 



In the case of hyperelastic materials, the stress-strain relation is furthermore assumed to derive from 
an elastic potential. Therefore the elasticity tensor is the second-derivative of the potential energy 
with respect to strain tensor. Thus, as a consequence of Schwarz's theorem, the hyperelasticity 
tensor possesses the major symmetry: 



The space of elasticity tensors is therefore the 21-D vector space Ela = S'^S'^(M'^). To each hypere- 
lastic material corresponds an elasticity tensor C but this association is not unique, there is a gauge 
group. Indeed, a designation based on the components of C in a fixed reference frame is relative to 
the choice of a fixed orientation of the material. If this problem is invisible for isotropic materials, 
it becomes more prominent as the anisotropy increases. For a fully anisotropic (triclinic) material 
how can it be decided whether two sets of components represent the same material ? To address 
this problem some approaches have been proposed in the literature so far. Roughly speaking, those 
propositions are based on the computation of the elasticity tensor spectral decomposition. But de- 
spite of being theoretically well-founded, the symmetry class identification relies on the evaluation 
of the roots multiplicity of the characteristic polynomial. Such a problem is known to be highly 
ill-posed [38] . Therefore this approach is difficult to handle in practice especially working with noise 
corrupted data. 

In order to avoid the computation of the spectral decomposition an alternative approach based 
on the invariants (and covariants) of the elasticity tensor is presently considered. Such a way to 
proceed is in the perspective of Boheler and al. [9] original work. It is interesting to notice that the 
general framework of the present paper is already well-known by the high energy physics community 



a' 







INVARIANT-BASED APPROACH TO SYMMETRY CLASS DETECTION 



3 



[1, 2, 33]. However, these methods do not seem to have been yet appHed in elasticity and especially 
for a representation as sophisticated as EI^(M^), which invariant algebra is not free. Reconstruction 
method based on this invariant-based approach, and numerical comparisons between approaches 
will be the objects of forthcoming papers. Before introducing in depth the proposed method, we 
wish to draw a short picture of this approach. Therefore, the next paragraph will be devoted to 
present the method in a nutshell. 

Invariant-based identification in a nutshell. As the material is rotated by an element g(= 
gj) G S0(3), the elasticity tensor C moves under the action of the 3D rotation group S0(3) on the 
space of elasticity tensors Ela 

Hence, from the point of view of linear elasticity, the classification of elastic materials can be 
assimilated to the description of the orbits of SO(3)-action on Ela. Practically, this is a difficult 
problem because the orbit space of a Lie group action is not a smooth manifold in general. Usually, 
the orbit space is a singular space which can described by the mean of the isotropic stratification [2] 
(described in section 1). In the field of elasticity, this isotropic stratification was first established 
by Forte and Vianello [19]^. 

Since the rotation group SO (3) is compact, the algebra of invariant polynomial functions on Ela 
is finitely generated [37] and separates the orbits [2, Appendix C]. Therefore, it is possible to find a 
finite set of polynomial invariants { Ji, . . . , Jn} which define a polynomial function J : Ela — >• 
such that 

j(c) = j(c') is a. 

This way, an SO(3)-orbit can be identified with a point in J(Ela) C M^. A classical result from 
algebraic geometry states that the image of a M-vector space under a polynomial map is a semi- 
algebraic set, that is a subset of M.^ defined by a boolean combination of polynomial equations 
and inequalities over M [13]. The orbit space is therefore a semi-algebraic set. The invariants 
{ Ji, . . . , Jn} are polynomials in 21 variables (dimEla = 21). To proceed towards our goal, we will 
make an extensive use of the decomposition of Ela into S0(3)-irreducible components, also known 
as the harmonic decomposition (c.f. section 3): 

Ela := S'^S^iR^) ~ 2EI°(M^) 2IH2(m3) IHI^(R^) 

where EI"(]R^) is the space of nth-order harmonic tensors on M^. Because of the lack of ambiguity, 
this notation will be shorten H". It is well known that dimH" = 2n -|- 1 for all n > 0. Historically 
this decomposition has already been used in the study of anisotropic elasticity tensors [6, 28, 7, 19]. 
We therefore refer to these references for a deeper insight into this topic. 

It is known [19] that Ela is divided into 8 conjugacy classes for the SO(3)-action. For each isotropy 
class [H], we build a linear slice Ela^ and we show that the equivalence relation C ~ C where 
C, C E ^[H] reduces to the linear action of a certain monodromy group on Ela^. The cardinal 
of is equal to the number of intersection points of a generic orbit in "^[h] with the linear slice 
Ela^. The closure of each stratum '^[h]j which we denote by ^[h] and refer to as a closed stratum, 
is a semialgebraic set which we describe explicitly by a finite number of equations and inequalities 
on the invariant polynomials { Ji, . . . , Jn}- For each class [H] for which F^ is finite, we give explicit 
relations between invariants to check whether C belongs to this class. 

Organization of the paper. The first section of this paper (section 1) is devoted to the introduc- 
tion of the framework uses for this study. There is no claim for any originality here, similar material 
can be found in [24, 33]. We have tried to make this paper accessible with only few prerequisites. 
The notions of stratification and linear slice are defined, together with a group theoretical method 
to compute its associated degree. In the second section (section 2), we illustrate these concepts on a 



The concept of symmetry class used by Forte and Vianello is equivalent to the notion of an isotropic stratum. 
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simple example: the stratification of the space of quadratic forms on under the SO(3)-action and 
we study the stratification of a n-uple of quadratic forms. This constitutes an extension of results on 
the stratification of 2 quadratic forms. Our results are provided in terms of harmonic tensors, which 
is not so usual. In section 3, the harmonic decomposition is introduced and the stratification of Ela 
discussed. The geometrical concepts introduced in section 1 will be enlighten in view of well-known 
results of anisotropic elasticity. In section 4, we study the strata of which have finite monodromy 
group. In each situation the stratification is established and algebraic criteria identifying symmetry 
classes are provided. These criteria are, up to authors' best knowledge, new. This section will be 
closed by a bifurcation diagram that sums up the symmetry breaking from H'^ together with the 
algebraic relations that induce those transitions. 

According to the restriction to linear slices, only 6 out of the 8 elastic symmetry classes pos- 
sess linear slices with finite monodromy group. Therefore, in the present paper, algebraic relations 
are only provided for the following anisotropic classes : isotropy ([S0(3)]), cubic ([O]), transverse 
isotropy ([0(2)]), tetragonal ([B4]), trigonal ([IDs]), orthotropic ([B2]). Nevertheless, the combi- 
nation of results from section 2 and section 4 completed by an assumption on the genericity of 
anisotropic of the elasticity tensor understudy (in a sense that will be defined in section 5.) allow 
the 8 symmetry classes to be identified. This paper will be concluded in section 6 and the extension 
of the method to a broader class of situations will be discussed. 

Main results. The main results of this paper are 

• In section 3 a rigorous geometric framework is provided that allows to recover important 
features of anisotropic elasticity [31]. Furthermore we draw a link between the conditions ap- 
pearing in the well-known Mehrabadi-Cowin theorem [16] and the elasticity tensor harmonic 
decomposition ; 

• In section 4, five sets of algebraic relations are provided to identify the following symme- 
try classes in H^: orthotropic ([B2]), trigonal ([B3]), tetragonal ([134]), transverse isotropic 
([S0(2)]) and cubic ([O]). Furthermore, we provide a parametrization of the corresponding 
strata by rational expressions involving up to 6 polynomial invariants, namely: J2, ■ ■ ■ , J7. 

• As an observation, we also prove that the invariants defined by the coefficients of the Betten 
polynomial [8] do not separate the orbits of Ela. 

Notations. To conclude this introduction let us specify some typographical conventions that will 
be used throughout this paper. Scalars will be noted with Greek letters ; minuscule letters indicate 
elements of H^, i.e deviator ; bold minuscule letters stand for elements of (8) M^, i.e. second order 
tensors ; elements of will be indicated by capital letters, meanwhile elasticity tensors (elements 
of ^^^^(M^)) win be noted by bold capital letters. Some minors exceptions to that rule may occur 
(such as for the stress a and strain e tensors) and will therefore be indicated in the text. 

1. Geometry of orbit spaces 
In this section, we consider a linear representation 

p : G^GL{V) 

of a compact real Lie group^ G on a finite dimensional R-vector space V. This action will be noted: 

g-v:= p{g)v 

Unless a real Lie group G acts without fixed point on V , the orbit space V/G is not a differentiahle 
manifold [2]. This is the main reason why the description of the orbit space is difficult in general. 
However, for a compact Lie group, the algebra of G-invariant polynomials, is known to separate the 



■In the following G will solely represent a compact real Lie group, therefore this precision will mostly be omitted. 
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orbit (c.f. subsection A. 5). Using this algebra, V/G can be described as a semialgebraic^ subset of 
M^. Let { Ji, J2, . . . , Jat} be a set of polynomial invariants which separate the G-orbits. Then, the 
map 

J -.v ^ {Ji{v), J2{v), ■■■ , Jn{v)) 

induces an algebraic homeomorphism between the orbit space V/G and J{V) C which IS a 
semialgebraic subset of M^. In particular, such a set is a stratified space. It is the union of manifolds 
of various dimensions, called strata., each of them being the union of orbits of a given isotropy type. 
This construction is the object of the present section. 

1.1. G-orbits. Two vectors vi and V2 are said to be G-related, and we write vi ~ V2 if there exists 
g G G such that f 2 = 9 ' vi. The set of all vectors v £ V which are related to v by G is called the 
G-orbit of v and is denoted by G ■ v: 

G-v:={g-v\geG} 
The G-orbits are compact submanifolds of V [2, 12]. 

1.2. Isotropy subgroups. The set of transformations of G which leave the vector v fixed forms 
a closed subgroup of G, which is called the isotropy subgroup (or symmetry group) of v. It will be 
denoted by G^: 

Gy := {g £ G; g ■ v = v} 
It can be shown that symmetry groups of G-related vectors are conjugate: 
Lemma 1.1. For any v £ V and any g £ G, we have 

Gg-v = g Gy g 

1.3. Isotropy classes. Two vectors vi and V2 are in the same isotropy class and we write vi ~ V2 
if their isotropy subgroups are conjugate in G, that is there exists g G G such that 

= 9 Gy^g 

Due to lemma 1.1, two vectors vi and V2 which are in the same G-orbit are in the same isotropy 
class but the converse is generally false. Therefore, the relation "to be in the same isotropy class" 
is weaker than the relation "to be in the same G-orbit". This observation is summed up by the 
following implication 

f 1 ~ W2 =^ I'l ~ V2 

The conjugacy class of a subgroup H will be denoted by [H]. The conjugacy class of an isotropy 
subgroup will be called an isotropy class. One can show that there is only a finite number of 
isotropy classes for a representation of a compact Lie group (see [12]). In the field of physics, for 
the SO(3)-action on tensor, this result is known as the Hermann theorem [26, 3]. On the set of 
conjugacy classes of closed subgroups of a G, there is a partial order induced by inclusion. It is 
defined as follows: 

[Hi] ^ [H2\ if Hi is conjugate to a subgroup of H2. 

Endowed with this partial order, the set of isotropy classes is a finite lattice^. In particular, it has 
a least element and a greatest element. 



semialgebraic set is a subset S of R" defined by a finite sequence of polynomial equations (of the form 
P{xi, . . . , Xn) — 0) and inequalities (of the form Q{xi, . . . , Xn)). 

^The lattice of conjugacy classes of all closed subgroups of SO(3) is described in Appendix B. 
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1.4. Isotropic stratification. The set of all vectors v in the same isotropy class defined by [H] is 
denoted '^[h] called a stratum. To avoid any misunderstanding it worths to note that a stratum 
is not a vector space, but a fiber bundle. There is only a finite number of (non empty) strata and 
each of them is a smooth submanifold of V [2, 12]. The partial order relation on conjugacy classes 
induces a (reverse) partial order relation on the strata 

The set of strata inherits therefore the structure of a finite lattice. An orbit G ■ v is said to be 
generic if it belongs to the least isotropy class. The generic stratum ^[Ho]j which is the union of 
generic orbits can be shown to be a dense and open set in V. Moreover, for each other stratum 
^[H]j we have dimSjj;^] < dimSfj^^jj]. The partition 

1/ = S[^„]US[^^]U---US[^„] 

is called the isotropic stratification of (y,p). '^[Ho] is the generic stratum. On the opposite side, the 
stratum which corresponds to the greatest isotropy subgroup is called the minimum stratum. 

Remark 1.2. Notice that the closure of a stratum ^[h]: denoted '^[h] ^^'i called the closed stratum 
associated to [H], corresponds to vectors v £ V such that v has at least isotropy [H], whereas ^[h] 
corresponds to vectors v £ V such that v has exactly isotropy [H]. 

1.5. Fixed point sets and normalizers. Let H be any subgroup of G. The set of vectors v € V 
which are fixed by H 

■.= {v£V\ h.v = V for all h £ H} , 

is called the fixed point set of H. For each v G , Gy C H. One can check that if Hi C H2 then 
(- yHi _ However, it may happen that V^'^ = V^^ but Hi 7^ H2. Furthermore, it should be 
pointed out that, in general, given Hi,H2 in the same conjugacy class [H], we have V^^ 7^ V^^. 
Notice also that a G-orbit belongs to an isotropy class greater than [H] iff it intersects the space 

Given a subgroup H of G, the normalizer of H, defined by 

N{H) ■.= {g£G\ gHg-' = H] 

is the maximal subgroup of G, in which H \s a normal subgroup. It can also be described as the 
symmetry group of H when G acts on its subgroups by conjugacy. 

Lemma 1.3. The space is N{H) -invariant. Moreover, if H = Gy^ is the isotropy group of 
some point vq G V, then N{H) is the maximal subgroup which leaves invariant . 

Proof. Let us first define 

K:= {g£G- g{V")=V^} 

Now let g G N{H) and v G . Given h £ H we have g~^hg G H and thus g~^hg - v = v. But then 
hg ■ V = g ■ V and since this is true for aW h £ H we get that g ■ v £ . That is is invariant 
under N{H), therefore N{H) C K. 

Suppose now that H = Gy^ for some point vq £ V . For allg £ K and h £ H,we have h-g-vo = g-VQ 
and hence g~^ ■ h ■ g ■ vq = vq. Therefore, g~^hg £ G„g = H. In other words, g~^Hg C H but, since 
this is also true for g~^, we have finally g~^Hg = H and therefore K C N{H). Hence K = N{H) 
which concludes the proof. □ 

1.6. Linear slices. One consequence of lemma 1.3 is that, the linear representation p : G — )• GL{V) 
induces a linear representation, of N{H) on , obtained by the restriction 

p^^H) : NiH) GL(F^). 



INVARIANT-BASED APPROACH TO SYMMETRY CLASS DETECTION 



7 



This induced representation is not faithful. However, when H is an isotropy group, its kernel is 
exactly H and we get a faithful linear representation 

PrH : — > GL{V^) where := N{H)/H. 

Furthermore, in that case, the equivalence relation defined by {V,p), when restricted to the set 
is the same as the equivalence relation defined by iV ■, Pn(h))- Indeed, let vi,V2 in 
, such that G^-^ = G^^ = H and V2 = g ■ vi. Then ghg~^ ■ V2 = V2, for each h ^ H. Thus 
ghg^^ £ Gy2= H and g G N{H). 

We have therefore reduced (locally, on the stratum ^[h]) the problem of describing the orbit 
space of V mod G by the orbit space of mod . This is especially meaningful when the group 
is finite, in which case, we say that is a linear slice. Then, each G-orbit intersect at 
most in a finite number of points and the natural map is a finite ramified covering 

of /r^ . The degree of is defined as the cardinal of F''^. It is equal to the index [N{H), H] 
of H in its normalizer N{H). It is also the number of points in a generic fiber of the covering 
— > /T^ . The group F^ is called the monodromy group of . 

1.7. Strata dimensions. The restriction of the projection map vr : 1/ — ?• V/G to a stratum "^[h] 

^/S[^j : ^ (1) 

is a fiber bundle with fiber G/H, i.e. a space that locally looks like Ti\^ij^/G x G/H. Therefore, we 
have 

dimS[H] = dim(G/F) + dim(S[^^]/G) 

where d\m[G/H) = dimG — dim//. 

Notice that the set n "^[h] is the subspace of vectors in which isotropy class is exactly 
[H]. It is an open and dense set in , i.e almost all elements of are in n '^[h]- These 
vectors correspond to the generic elements in . 

Now, since the equivalence relation defined by {V,p), when restricted to the space , is the same 
as the equivalence relation defined by {V^ , p-ps), the basis H^h^^/G of the bundle (1) is diffeomorphic 
to {y^ n S[j:^])/F^ (see [12] for a more rigorous justification). Thus 

d\^T.^H]/G = dim(y^ n S[//])/F^ = dhuV^ - dimF^ 

dimSfj:^] = dimy^ + dim(G/iV(F)) = dimT/^ + dimG - dimiV(F) 

Mechanical illustrations of the physical contents of those geometric concepts will be given in 
subsection 3.2. 

Remark 1.4. The dimension of can be computed using the trace formula for the Reynolds 
operator (see [24]). Letting xp be the character of the representation iy,p), we have 

di^V^ = ±-S2xp{h) (3) 

if // is a finite group (and the preceding formula has to be replaced by the Haar integral over H for 
an infinite compact group). Explicit analytical formulas based on (3) were obtained in [4] and [24]. 
They can be found in Appendix C. 



ramified covering is a generalization of the concept of a covering map. It was introduced in the theory of 
Riemann surfaces and corresponds to the situation where a finite (or discrete) group acts on a manifold but with fixed 
points. A typical example is given by the map 2 i-^ z" on the complex plane. Outside 0, this map is a covering but 
is a singularity, a point of ramification. In higher dimension, the set of ramifications points may be more complicated. 
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1.8. Implicit equations for closed strata. Let [H] be a fixed isotropy class and ^[h] C V the 
corresponding stratum. Given a finite set of invariant polynomials which separate the orbits of G , 
the closed stratum '^[h] can be characterized by a finite set of relations (equations and inequalities) 
on { Ji, J2, . . . , Jat}. To obtain these equations, we fix a subgroup H in the conjugacy class [H] 
and choose linear coordinates {x^)i<i<q on . Then, we evaluate {Ji, J2, ■ ■ ■ , Jn} on (as 
polynomials in the x*) and we try to obtain implicit equations on the for the parametric system 

Jk = Pk{x\...x''), k = l,...,N, PkeR[Xi,...,Xg] (4) 

satisfied by the restriction of { Ji, J2, . . . , Ja?} to . 

This is a difficult task in general (one might consider [17] for a full discussion about the im- 
plicitization problem). Moreover, it could happen that the algebraic variety Vj defined by such an 
implicit system is bigger than the variety Vp, defined by the parametric system (4). Besides, we 
are confronted to the difficulty that we work over M which is not algebraically closed. To overcome 
these difficulties, we observe that the restrictions of { Ji, J2, . . . , Jat} to are F^-invariant and 
can therefore be expressed as polynomial expressions of some generators fii , . . . , o",. of the invariant 
algebra of the monodromy group on . This observation reduces the problem to consider first 
the implicitization problem for the parametric system 

Jk = PkicTi,. . . ,ar), k = l,...,N, pk £R[Xi, . . . ,Xr] (5) 

To solve this problem, we use a Groebner basis [17] (see subsection A. 4) or a regular chain[ll]. The 
main observation is that, in each case we consider, the system (5) leads to a system of relations 
(syzygies) which characterizes the closed stratum '^[h] 

Sji,h,J2,...,JN) = 0, j = l,...l Sj £R[Xi,...,Xn], (6) 

and a system 

ai = Ri{Ji,J2,...,JN), i = l,...,r Ri £ R{Xi, . . . , Xn), (7) 
which express the crj as rational^ functions of { Ji, J2, . . . , Jn} on the open stratum ^[h]- Beware, 
however, that these rational expressions may not be unique. However if some can be written as 
Pi/Qi as well as P2/Q2 then P1Q2 — P2Q1 belongs to the ideal generated by the Sj. 

Because the solutions are rational, the fact that the field on which we work is real or complex 
does not matter at this level. Therefore, for each real solution (Ji, J2, . . . , Jn) of (6), it corresponds 
a unique real solution (cri,... ,ar) given by (7). Nevertheless, we cannot take for granted that a 
real solution (cti, . . . , ar) of (7) corresponds to a real point {x^, . . . , x'^) in . For this, we need to 
compute an additional system of inequalities on the Uj, or equivalently on the Jk, which permits to 
exclude complex solutions of 

ai = Qi{x\...x'i), i = l,...,r, Qi e R[Xi, . . . , Xq] 

However, in the present paper for the linear slices we consider, the only monodromy group that we 
encounter are^ 1, S2 and 63 and the action is the standard one in appropriate coordinates. In these 
cases, and more generally when the monodromy group is isomorphic to S„ and acts by permutation 
on the coordinates, the invariants fii, . . . , (T„ are algebraically independent, r = q = n and x^, . . . 
are the roots of the polynomial 

Therefore, the problem reduces to find conditions on the ai that ensure that all the roots of p are 
real. The solution is due to Hermite, [13]. He has proved that the number of distinct real roots of a 
real polynomial p of degree n is equal to the signature of the Hankel matrix B{p) := (S'j+j_2)i<j 



The fact that the solutions are rational will not be justified here. We just observe that this is the case for all 
classes we have treated so far in this article. 

'''where &„ is the the group of permutations acting on n elements (symmetric group). 
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where Sj. := 'Y17=ii-'^^)^ power sum of the roots of p. In particular, p has real roots if and only 

if B(p) is non- negative. 

In the forthcoming section this method will be applied to establish the stratification of a n-uple 
of quadratic forms. Besides being just an illustration, the results obtained will be of great interest 
for the symmetry detection of elasticity tensors and used in section 3 and section 5. 

2. The isotropic stratification of a n-UPLE OF QUADRATIC FORMS 
In this section, we will study the isotropic stratification of the space 

n 
k=l 

For n = 2, that is for a couple of quadratic forms, we get the lower order terms of the harmonic 
decomposition of Ela (c.f. section 3). Invariant conditions which characterize each symmetry class 
will be given. 

2.1. The isotropic stratification of one quadratic form. The lattice of isotropy classes for the 
tensorial representation of S0(3) on the space of symmetric 3-matrices, /^^(M^^^ jg ^ total order 

[B2]=<[0(2)]^[SO(3)]. 

Since the least isotropy class, [B2], is not trivial, the fixed point set S'2(M3)ID'2 (which corresponds 
to the subspace of diagonal matrices) is a global linear slice. Each SO(3)-orbit intersects this 
slice (a symmetric matrix is diagonalizable in a direct orthogonal frame). This linear slice is of 
dimension 3 and degree 6 (see Appendix B). Generically, an orbit intersects this slice in exactly 6 
points which correspond to the permutation of the eigenvalues Ai , A2 , A3 . The monodromy group 
N{l\}2)/02 ^ &3 acts on {Ai,A2,A3} by permutation. The stratification problem is reduced to 
the computation of the number of distinct eigenvalues, namely (3,2,1). The invariant algebra of 
(S'2(R3)''^2^ ig ^j^g fj-QQ algebra generated by the 3 symmetric elementary functions o"i,(T2,o"3. 
According to Hermite (see [13]) the characteristic polynomial of a has 3 distinct real roots iff the 
quadratic form 

/ 3 2ai (72 \ 

5(a) := 2ai 2ai^-2a2 (Jia2-2,(Tz 

\ (T2 (71(72-3(73 - 2 (Ji(73 + (72^ / 

is positive definite. The three principal minors of B are given by 

' Ai = 3, 
< A2 = 2al - 6fJ2, 

^ A3 = -27(73^+ (18(710-2 -4 (71^) (73 - 4 (72^ + (71^(722. 

These minors are precisely those which have been introduced by Hermite to count the number of 
real roots of a real polynomial of degree 3. A3 is the discriminant of the polynomial 

A3 = (Al-A2)'(A2-A3)'(Al-A3)^ 

whereas 

A2 = (Ai - A2)' + (A2 - A3)2 + (Ai - A3)' 

vanishes when all the roots are equal. They give directly the classification we are looking for, which 
is summarized in table 1. 
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H 


N{H) 




card r-'^ 


strata T^h 


dim V" 


dimS[j|/j/G 


dim TiH 


D2 




63 


6 


A3 >0 


3 


3 


6 


0(2) 


0(2) 


1 


1 


A3 = 0, A2 > 


2 


2 


4 


S0(3) 


S0(3) 


1 


1 


A3 = 0, A2 = 


1 


1 


1 



Table 1. Isotropy classes for the representation of S0(3) on symmetric 3-matrices. 



H 


N{H) 




card 


dim 


dim 


dimS[^] 


1 


S0(3) 


S0(3) 


00 


6n 


6n — 3 


6n 


Z2 


0(2) 


0(2) 


00 


in 


4n - 1 


4n + 2 


B2 





63 


6 


3n 


3n 


3n + 3 


0(2) 


0(2) 


1 


1 


2n 


2n 


2n + 2 


S0(3) 


S0(3) 


1 


1 


n 


n 


n 



Table 2. Isotropy classes for the representation of S0(3) on a n-uple of quadratic forms. 



2.2. The isotropic stratification of n quadratic forms. In this subsection, we consider a n- 
uple of quadratic forms (ai, . . . ,a„). The lattice of isotropy classes is the following totally ordered 
set (see [32] for a general algorithm to compute isotropy classes) 

[1] < [Z2] < [B2] < [0(2)] < [S0(3)] 

The stratification is summarized in table 2. To establish the bifurcation conditions we need to 
introduce the commutator of two quadratic forms. Given two quadratic forms (a, b) on the euclidean 
space M^, the commutator 

[a, b] := ab — ba 

is a skew-symmetric endomorphism $7, which can be represented, in the oriented vector space M^, 
by a vector oj. We will denote by a; 1— )• := this linear isomorphism between vectors of M'^ and 
skew-symmetric second order tensors on M^. Recall that for every g G S0(3) and to € M^, we have 

j(gw) = gj(^)g*- 

In particular 

.;(a,b) :=r'([a,b]) 
is a covariant of degree 2 of the pair (a, b). 

Theorem 2.1. The characterization of isotropy classes in @^=qS'^{E?) is summarized below: 

(1) A n-uple (ai, . . . ,a.„) has at least isotropy [Z2] iff there exists a non-vanishing vector cu € 
and real numbers Xk,i such that 

w(ak,ai) = Afc,;tJ, l<k,l<n, 
akW Aa; = 0, 1 < k < n, 

where A is the classical vector product ofW^. 

(2) A n-uple (ai, . . . ,a„) has at least isotropy [D2] iff 

a;(ak,ai) = 0, l<k,l<n. (9) 

(3) A n-uple (ai, . . . ,a„) has at least isotropy [0(2)] iff 

25[tr(4)]' = 54[tr(4)]^ 1 < k < n, ^^^^ 

[tiiakai)f = [tr(4)] [tr(af)] , 1 < fc, / < n. 



INVARIANT-BASED APPROACH TO SYMMETRY CLASS DETECTION 



11 



where = dev(afc) is the deviatoric part o/a^. 
(4) A n-uple (ai, . . . ,a„) has isotropy [S0(3)] iff 

ak = 0, l<k<n, (11) 

where = dev(afc) is the deviatoric part o/a^. 

Proof. (1) If a n-uple (ai, . . . ,a„) has at least isotropy [^2], its elements share a common principal 
axis (the axis of the non trivial rotation a in Z2). Let a; be a unit vector parallel to this axis. Then 
uj belongs to the kernel of the commutator of each pair (ak,ai) and thus a;(ak,ai) is proportional 
to 00. Relations (8) are a consequence of these facts. Conversely, if relations (8) are satisfied, then 
w is a common eigenvector of the n-uple (ai, . . . ,a„). In that case, the isotropy group of the pair 
(a, b) is at least in the class [Z2]. 

(2) If a n-uple (ai, . . . ,a„) has at least isotropy [O2], then there exists a g £ S0(3) which brings 
all the afc to a diagonal form. Thus they commute. Conversely, if the a^ commute, they can be 
diagonalized simultaneously and the isotropy group of the n-uple (ai, . . . , a„) is at least in the class 

m. 

(3) A n-uple (ai,...,a„) has at least isotropy [0(2)] iff all the a^ can be diagonalized simul- 
taneously and, moreover, iff their respective deviator has (in a well chosen basis) a diagonal form 
proportional to diag(l, 1, —2). The first relation of equations (10) expresses that the deviator of a^ 
has a double eigenvalue, the second relation expresses that the deviators and ai are proportional 
(Cauchy-Schwarz relation). Conversely, if equations (10) are satisfied, then all the deviators 
commute and each of them has a double eigenvalue. It is then possible to find a basis in which 
all the ttfc are proportional to diag(l, 1, —2). Therefore the n-uple (ai, . . . ,a„) has at least isotropy 
[0(2)]. 

(4) The last assertion is trivial. □ 

3. The space of elasticity tensors 

3.1. The harmonic decomposition of elasticity tensors. To establish the isotropic stratifica- 
tion of a tensor space, the first step is to decompose it into a direct sum of elementary pieces. For the 
SO(3)-action, this decomposition is known as the harmonic decomposition. The S0(3)-irreducible 
pieces are symmetric and traceless tensors (i.e. harmonic tensors). The space of elasticity tensors 
admits the following harmonic decomposition [6, 7, 19, 21]: 

Ela ~ 211° e 2M^ e (12) 

Therefore, each C € Ela can be written as C = (A, /U, a, 6, D) where X, fi € M^, a, 6 G and G 
and such that for all g G SO (3) 

C = p{g) • C A = A, /X = /X, a = p2{g) ■a,b = p2{g) ■ b, D = pi{g) ■ D 

Elements of the set {X, p,a,b, D) are covariants (c.f. subsection A. 2) to C : (A,/i) are invariants, 
(a, b) are /92-covariants, and D is p4-covariant. The explicit harmonic decomposition is well known, 
given an orthonormal frame (ei, 62, 63) of the Euclidean space, we get [9]: 

Cijki = X6ij6ki + p{SikSji + dii6jk) + SijOki + Shmj + 6ikbji + 6jibik + Subjk + Sjkbu + AjH (13) 

where the metric is written qij := 5ij. This formula can be inverted to obtain the 5 harmonic tensors 
(A, fi, a, b, D) from C. On Ela = ^^^^(M^), there are only two different traces: 

3 3 
dij = (tri2 C)ij := ^ Ckkij, Vij = (tri3 C)ij := ^ Ckikj- 
k=l k=l 

where d and v are known as, respectively, the dilatation tensor and the Voigt tensor [16, 14]. 
Starting with (13) we get 

d = (3A + 2/i)q + 3a + 46, v = (A + 4^)q + 2a + 56. (14) 
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Taking the traces of each equation, one obtains: 

tr(d) = 9A + 6/u, tr(v) = 3A + 12/i 

and, finahy: 

A = ^(2tr(d) - tr(v)), fi = ^(- tr(d) + 3tr(v)) 

a = i(5dev(d) -4dev(v)), b= ^(-2 dev(d) + 3 dev(v)) 

where we write dev(a) := a — | tr(a) q the deviatoric part of a 2nd-order tensor. 

The isotropy classes of elasticity tensors have been computed in [19, 10]. There are exactly 8 
classes: the isotropic class [S0(3)], the cubic class [O], the transversely isotropic class [0(2)], the 
trigonal class [D3], the tetragonal class [D4], the orthotropic class [B2], the monoclinic class [Z2] 
and the totally anisotropic class [1] (also know as the triclinic class). We refer to Appendix B for 
the precise definitions of these groups. This stratification is a consequence of the index symmetries 
of elasticity tensors. It is not a general result that encompass all kind of 4th-order tensors. For 
example, the photo-elasticity tensors space which is a 4th-order tensor space solely endowed with 
minor symmetries is divided into 12 symmetry classes [20]. 

3.2. The isotropic stratification of Ela. The stratification of Ela is provided in the table 3. The 



H 


N{H) 




card r-'^ 


dim 


dim Yi^fj^/G 


dimS[^] 


1 


S0(3) 


S0(3) 


CO 


21 


18 


21 


Z2 


0(2) 


0(2) 


00 


13 


12 


15 


©2 





63 


6 


9 


9 


12 


D3 


Be 


62 


2 


6 


6 


9 


D4 


Dg 


©2 


2 


6 


6 


9 


0(2) 


0(2) 


1 


1 


5 


5 


7 








1 


1 


3 


3 


6 


S0(3) 


S0(3) 


1 


1 


2 


2 


2 



Table 3. Isotropy classes for Ela. 



dimensions of ^[h] ^-iid are obtained using the formulas (2) of subsection 1.7. And dimV^ 

was computed using the trace formula (3) and the explicit formulas provided Appendix C. Table 3 
sums up the different dimensions that are classically associated to anisotropic elasticity tensors in 
the literature. All these dimensions are linked by the formula 

dim S = dim S [j/] /G + dim + en, 

where Ti^fj^/G is the orbit space for tensors in isotropy class [H], = N{H)/H is the monodromy 
group and en = A\m.G/N[H). Notice that en — [H] is the set of distinct subgroups of G which are 
conjugate to H. These results of elasticity enlighten the physical contents of the geometric concepts 
introduced in section 1. 

: The fixed-point set is the vector space of tensors having at least H as isotropy group. 
Its dimension equals the dimension of the matrix space used to represent them. The space 
can be constructed for any closed subgroup H, not only for an isotropy subgroup. In 
elasticity, for instance, fixed-point sets have been constructed for Z3, Z4 (which are not 
isotropy subgroups). This is why it was believed for a long time that there exists 10 type 
of elasticity. The construction of such matrix spaces is detailed, for example, in [15] for 
classical elasticity, and extended to second order elasticity in [5] . It worthwhile emphasizing 
the fact that for a non-isotropy subgroup, is solely a fixed-point set and in no way a 
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slice. It is geometrically meaningless for the description of the stratification, that is the 
geometry of the orbit space. 

'^[H] '• The set of tensors which symmetry group is conjugate to H forms a stratum, or a 
symmetry class according to Forte and Vianello [19]. We got the following dimension relation 
dimE[j|^] = dimV^ + eH where en = dimS0(3)/A^(-fr) is the number of Euler angles needed 
to define the appropriate coordinate system in which the elastic tensor has the minimum 
number of elastic constants [31]. The strata dimension corresponds to what Norris, for 
example, calls the number of independent parameters characterizing a tensor. 

T,[jj^/G : This is the orbit space of '^[h]i i-^- each point of this space represents the G-orbit 
of a tensor which symmetry group is conjugate to H. The intrinsic properties of a tensor 
are characterized by its orbit. In a slice the image of the orbit is given by monodromy 
group, and we have the dimensional relation diuiV^ = dimSjji^j/G + dimF^. Therefore 
when is discrete, Tj^fj^/G and have the same dimension, this is the situation for most 
of elasticity symmetry classes. At the opposite, when is continuous, the dimension of 
the orbit is strictly smaller than the one of the slice, in elasticity this case occurs for 1,^2. 

Remark 3.1. It is known for long time that for the symmetry classes 1,^2, there is no canonical 
way to obtain a basis in which the number of non-zero components is minimal. More precisely, an 
element of the orbit of such a tensor in an appropriate fixed point set does not have necessarily 
the minimal number of non-zero components and thus, such a basis has to be computed arbitrarily. 
This is due to the fact that in these cases, is continuous. For the trivial class, [1], ~ V and 
there is no natural way to obtain a normal form. The monodromy group F^ ~ SO (3) and there are 
three arbitrary degree of freedom to choose a nice basis. For the other cases, is smaller than 
V but the number of non-zero components of a tensor in can usually be lowered arbitrarily. 
The monodromy F^ ~ 0(2) and there is one arbitrary degree of freedom to lower the number of 
non-zero components. In each cases, it appears clearly that one can properly choose a basis in 
which those angular parameters equal zeros, such a basis is non necessarily unique. In such a basis 
the number of independent parameters is lowered and we obtain the minimum number of elastic 
constants, as discussed, for example, by Norris [31]. This minimum number of elastic constants is 
equal to the orbit dimension. 

The results appearing in table 3 are well-known in elasticity, but their constructions are usually 
ad-hoc and do not stand on rigorous mathematical foundations. Therefore, we hope to provide such 
a firm geometric construction. Conversely, as the geometrical concepts we use are abstract, these 
well-known mechanical results help to illustrate their physical contents. As it has been said before 
the space of elasticity tensors can be decomposed in the following way: 

Ela ~ 2]HI° e 2lf e (15) 

The section 4 will be devoted to the stratification of H^. The stratification of 2EI^ © 2EI^ is a 
consequence of the results obtained in the subsection 2.2 in the case n = 2. In such a case a direct 
link between the bifurcation conditions given (2.1) and the ones appearing in the Cowin-Mehrabadi 
theorem can be drawn. 

3.3. The Cowin-Mehrabadi condition. The information contained in the Voigt and dilatational 
tensors (14) solely concerns the covariants (A,//, a, 6). As a consequence any statement on d and v 
can be translated in terms of (A, /i, a, b), and conversely. 

The condition (8) for (a, b) to have at least isotropy [Z2] means that uj = [a, b] is a common 
eigenvector of (a, b). If we note n = w, it is a straightforward computation to show that this 
condition on (a, b) is equivalent to 

dijUj = 6ni, and VijUj = urii. 
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These conditions are nothing but the two first ones of the Cowin-Mehrabadi (CM) theorem for the 
existence of a symmetry plane. Let summarize this theorem as it is formulated in [36]. 

Theorem 3.2 (Cowin-Mehrabadi). Let n he a unit vector normal to a plane of material symmetry, 
if m is any vector on the symmetry plane a set of necessary and sufficient conditions for n to be 
normal to the symmetry plane is 

Cijkknj = {Cpqttnpnq)ni, 

Cisks'^k — {Cpqrqripny^ni, 
Cijks'^j'^s'^k — {CpqptTlpTlqTlj^Tli^Tli^ 
Cijksf^jlT^s'lT'k — iCpqrtn'pn^qn'r^t)^i- 

The two first conditions of this theorem can be rephrased in terms of the second order covariants 
(a, b) of the harmonic decomposition. The bifurcation condition (8) provides a method to find a 
good candidate to be the normal of a symmetry plane uj = [a, b]. In the same spirit the third and 
fourth conditions are indeed conditions on D. Therefore the theorem can be reformulated in term 
of harmonic covariant using the condition (8): 

Theorem 3.3 (Cowin-Mehrabadi). Let C = {X, fi,a,b, D) be an elasticity tensor. Let suppose there 
exist a non null-vector uj = [a, b], this vector is a vector normal to a plane of material symmetry, if 
there exists uj-^ verifying oj-^ ■ cj-*" = and such that the following conditions are satisfied 

UJ f\{D : 00 ® uj)(jj = 0, 

UJ h{D : uj^ ®uj^)uj = 0. 

4. The isotropic stratification of H"^ 

In this section, the isotropic stratification of will be established. This space is isomorphic to the 
higher order component of the harmonic decomposition of Ela. Invariant conditions characterizing 
each finite monodromy symmetry class will be given. These results provided new sets of necessary 
conditions to identify higher order symmetry classes, and will be used in section 5. 

4.1. Lattice of isotropy. The isotropy classes of the 9-D space are the same as the 8 classes 
of Ela. The corresponding lattice® is illustrated in figure 1. With the convention that a group at 
the starting point of an arrow is conjugate to a subgroup of the group pointed by the arrow. 

Contrary to what happens with quadratic forms^, there is no nice description of the isotropic 
stratification. The aim of this section is to describe the strata which have finite monodromy group 
r-'^: the cubic class [O], the transversely isotropic class [0(2)], the trigonal class [D3], the tetragonal 
class [P4] and the orthotropic class [B2]. For each of them, explicit results will be provided. The 
algebraic procedure used to obtain their solutions was explained in subsection 1.8. 



It should be noticed that it is not possible to realize this lattice of conjugacy classes as a lattice of subgroups (for 
the inclusion relation). In particular, if each conjugacy class in figure 1 is replaced by the corresponding subgroup 
defined in Appendix B, then almost all the arrows corresponds to an inclusion relation but not the arrow from (4) to 
(2) : D3 is not a subgroup of O but is conjugate to a subgroup of O. 

^It is quite remarkable that for quadratic forms in every dimension n, the generic stratum has non trivial isotropy 
and that the corresponding linear slice has finite monodromy 6„. This is the main reason why every orbit has a 
natural normal form, namely a diagonal representation (up to the finite monodromy group In comparison, the 

main difficulty in linear elasticity is the fact that for a generic elasticity tensor, the isotropy group is trivial and there 
is no natural normal form. 
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@ isotropic: [S0(3)] 



transversely ii 




(S) triclinic: [1] 

Figure 1. The lattice of isotropy classes for W^. 



H 


N{H) 




card 


dim 


dim T,^jjj/G 


dimS[^] 


1 


S0(3) 


S0(3) 


00 


9 


6 


9 


Z2 


0(2) 


0(2) 


00 


5 


4 


7 


©2 





63 


6 


3 


3 


6 


D3 


Be 


62 


2 


2 


2 


5 


D4 


Bg 


62 


2 


2 


2 


5 


0(2) 


0(2) 


1 


1 


1 


1 


3 








1 


1 


1 


1 


4 


S0(3) 


S0(3) 


1 


1 












Table 4. Isotropy classes for the tensorial representation of S0(3) on H^. 



4.2. Global parametrization. On an euclidean vector space, any bilinear form can be represented 
by a symmetric linear operator. This applies in particular to an elasticity tensor C which may 
be considered as a quadratic form on S'^(M^), where the Euclidean structure is induced by the 
standard inner product on M^. It is therefore customary to represent an elasticity tensor C by its 
corresponding symmetric linear operator C on S''^(M^) (relatively to the canonical scalar product 
(a, 6) = tr(a6)). The corresponding matrix representation[35, 15] is given by 



C = 



/ Cii 


C12 


Cl3 


\/2ci4 


V2C15 


\^Ci6 


\ 


C12 


C22 


C23 


V2C24 


V2C25 


V2C26 




Cl3 


C23 


C33 


\/2C34 


\/2C35 


\/2C36 




\/2ci4 


\/2C24 


\/2C34 


2 C44 


2 C45 


2 C46 




V2C15 


V2C25 


V2C35 


2 C45 


2 C55 


2 C56 




V V2C16 


\/2C26 


\/2c36 


2 C46 


2 C56 


2 C66 


/ 



where Cmn are the components of the elasticity tensor in an orthonormal frame. We use the conven- 
tional rule to recode a pair of indices (i, j) = 1, 2, 3) by a integer m = a{i,j) (m = 1, 2, . . . 6), 
where a{i,i) := z for 1 < i < 3 and a{i,j) := 9 — (i + j) for z / j. 
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Given a general tensor D S H^, we denote by D. the corresponding symmetric linear operator. In 
the canonical basis of M^, D_ has the following matrix representation, depending on 9 real parameters: 



D 



V2DT 



12 



2D22 



where ^ represents the transposition, 



D 



11 




hg 
-hg - h-j 



and 



D 



22 



hs 

-/is - hy 



-hi - /13 

-/l2 - /14 



D 



12 




h2 
-h2 - hi 

/l4 



hi 

-hi - hs. 



-hi - hs 




-/15 - he 

Using this parametrization, D has the following invariance properties: 

(i) Z2-invariant iif /i2 = /14 = /is = /le = , 

(ii) D2-invariant iif /i^ = pour 1 < k < 6 , 

(iii) Da-invariant iif /ii = /i2 = /13 = /i4 = /le = et /17 = /ig = —4:hg , 

(iv) D4-invariant iif h^ = pour 1 < /c < 6 et /17 = /ig , 

(v) 0(2)-invariant iff /i^ = pour 1 < /c < 6 et /17 = /ig = — 4/19 , 

(vi) O-invariant iff /i^ = pour 1 < /c < 6 et /17 = /ig = /19 , 

(vii) S0(3)-invariant iff /i^ = pour 1 < k < 9 . 

In the following subsections, the fundamental invariants of 4-th order harmonic tensors will be 
used. As these results are not well-known we will, in the next lines, summarize their expressions 
as obtained by Boelher and al.[9]. For technical details and historical considerations, we therefore 
refer to this publication, and the references therein. 



Proposition 4.1. Let D G 



d2 
dg 



The 2nd-order tensors d2, . . . , diQ.' 
2^ - tri3(Z)3) 



triap^) dg 
d2Dd2 de 
d|L>2d2 d 



9 



d3 



d4 = 
d7 = 

dio 



di 



d^Dd2 



(16) 



are covariant to D. The nine fundamental invariants are defined by: 

Jfc:=tr(dfc), A; = 2,..., 10. 

The first 6 invariants J2,..., J7 are algebraically independent. The last 3 ones Jg, Jg, Jio are 
linked to the formers by polynomial relations. These fundamental syzygies were computed in [34]. 

Remark 4.2. Notice that the first invariant J2{D) = {D, D) is the squared norm of D, it corresponds 
to the squared Frobenius norm of D- In particular, J2{T)) = if and only if D = 0. 



4.3. Cubic symmetry 

tion: 



D 



A O-invariant tensor D G has the following matrix representa- 



(17) 





85 


-A5 


-4 6 















-A5 


86 


-4 6 















-A5 


-4 6 


86 
























-86 
























-86 







1 

















-86 


) 
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where 5 € M. The evaluation of the invariants (16) on this slice gives the following parametric 
system: 



J2 


= 480<5^ 


Js 


= 1920 5^ 


J4 


= 76800 5^ 




= 0, 


Je 


= 12288000 5^ 


Ji 


= 0, 




= 0, 


J9 


= 0, 




= 0. 



We notice that the parameter 4(5 = J3/J2 is a rational invariant. The linear slice corresponding to 
the octahedral group is of degree 1, which was already known from the algebraic argument explained 
in Appendix B. 

Proposition 4.3. A harmonic tensor D G is in the closed stratum SfQ] if and only if the 
invariants J2{D) ■ ■ ■ Jio{D) satisfy the following polynomial relations 

3J4 = J2', J5 = 0, 30J32 = J2^ 9J6 = J2^ 
J7 = 0, J8 = 0, Jg = 0, Jl0 = 0. 

D is in the cubic class [O] if moreover J2{D) 7^ 0. In that case, it admits the normal form (17) 
where 45 := Jz{D)/J2{D). 

We observe that VJ2, J3 G M, 5 G M. Therefore we do not need to add any inequality to the 
syzygy system (18). 

Remark 4.4. All 2nd-order tensors in Wi? covariant to D vanish. 

4.4. Transversely isotropic symmetry ([0(2)]). A 0(2)-invariant tensor G has the fol- 
lowing matrix representation: 



D 



/ 


35 


5 


-45 











\ 




5 


35 


-45 















-A5 


-45 


85 
























-8 5 
























-85 







V 

















25 


/ 



(19) 



where 5 G M. The evaluation of the invariants (16) on this slice gives: 

J2 = 280 5^ Ja = 720 5^ J4 = 32800 5^ , 

J5 = 80000 5^ Je = 4528000 5^ J7 = 17600000 5^ , 

Jg = 211200000 5^ Jg = 3872000000 5^ , Jio = 46464000000 5^° 

The parameter 5 = 7 J3/I8 J2 is a rational invariant. The linear slice corresponding to the group 
0(2) is of degree 1, which was already known (see Appendix B). 

Proposition 4.5. A harmonic tensor D G H'^ is in the closed stratum S[o(2)] ^/ only if the 
invariants J2{D) " " " '^io(^) satisfy the following polynomial relations 

98 J4 = 41 J2^ 63 J5 = 25 J3 J2, 3430 Js^ = 81 J2^ 1372 Jg = 283 J2^ 

9 4 ^ (20) 

882 J7 = 275 J2 V3, 4802 Jg = 165 J2^, 12348 Jg = 3025 J2 V3, 67228 Jio = 1815 J2^ 

D is in the transversely isotropic class [0(2)] if moreover J2{D) 7^ 0. In that case, it admits the 
normal form (19) where 5 = 7 J^{D) / IS J2{D) . 

As for [O] the expression of 5 ensures that, for all J2, J3 G M, the solution 5 is real. Therefore no 
inequality is needed. This fact is indeed general to any linear slice of degree 1. 

Remark 4.6. All 2nd-order tensors in 'M? covariant to D are multiple of diag(l, 1, —2) and are 
eigenvectors of D corresponding to the eigenvalue 125. 
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4.5. Trigonal symmetry ([B3]). A IDs-invariant tensor D G has the following matrix repre- 
sentation: 



D 



/ 


35 


S 


-45 


-V2a 








\ 




6 


35 


-45 


V2a 












-45 


-4 5 


85 

















V2a 





-8 5 
























-85 


-2a 




V 














-2a 


25 


/ 



where {6, a) G M?. The action of the monodromy group ©2 is given by: 

5 I—)- 5, o" I— 7- —a . 

The evaluation of the invariants (16) on this slice gives: 

J2 = 280 5^ + 16 

J3 = 144 5 (552 -cj2) 

J4 = 32800 5*^ + 2720 a'^6'^ + 88 cj^ 

J5 = 32 5 {-a^ + 50 6^f 

Jg = 4528000 5^ + 436800 £7^5^ 20640 0-^5^ + 496 ct^ 
J7 = 320 5 (22 5^ + cj2) (-ci^ + 50 5^)^ 
Js = 3840 52 (22 5^ + (t2) (-a^ + 50 5^)^ 



'8 

J9 



3200 5 f-a^ + 50 52)' (2252 +ct2)' 



2\2 



38400 52 (-a^ + 50 5^) (22 5^ +cj2) 
We notice that the parameter: 

5 



1 



J5 



(21) 



4 Js^ - 3 J4 

is a rational invariant, therefore for all J2,J4,J5 G M, the parameter 5 is real. However, the minimal 
equation satisfied by a is of degree 2 



Jo = 280 52 + 16 . 



(22) 



A condition is needed for a to be real. This condition, which is > 0, is given by the following 
inequality on J2 , J4 , J5 : 



2J2(J| - 3 J4)2 - 35 J| > 



(23) 



Having computed 5, the equation (22) has two roots with opposite signs. The linear slice corre- 
sponding to the group D3 is of degree 2, which was already known (see Appendix B). 
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Proposition 4.7. A harmonic tensor D G is in the closed stratum Sp^] if and only if the 
invariants J2{D) ■ ■ ■ Jiq{D) satisfy the following polynomial relations 



192 Je 
36 J7 
768 Ji^ 
240 Js 
576 J4 J5 
1152 Jg 
1440 Js^ 
8640 Jio 



-51J2'' + 216J2J4 + 10J3^ 

-2 J2V3 + 6J3J4 + 27J2J5 
-99 J2^ + 552 J2 V4 + 10 J2 Js^ + 240 J3 J5 
-33 J2^ + 96 J2 ^4 + 30 J2 Js^ + 40 J3 J5 
-41 J2 V3 + 120 J2 J3 J4 + 216 J2 V5 + 30 J3^ 
-99 J2 V3 + 296 J2 J3 J4 + 648 J2 V5 + 10 J3^ 



(24) 



-11 J2^ 32 J2 ^4 - 70 72^3^ + 240 J2 J3 J5 + 240 J3 V4 

-891 J2^ 2592 J2 V4 + 730 J2 ^3^ + 2160 J2 J3 J5 + 240 J32 J4 

together with inequality (23). D is in the trigonal class [B3] if moreover 3 J4 - J| / and 98 J4 - 
41J| / 0. In that case, it admits the normal form (21) where 

= — T ^5 

4 J22 - 3 J4 

and where a is the positive root of J2 = 280 5^ -|- 16 a^. 

Remark 4.8. All 2nd-order tensors in covariant to D are multiple of diag(l,l,— 2) and are 
eigenvectors of D corresponding to the eigenvalue 125. 

1)4]). A D4-invariant tensor D £ has the following matrix rep- 



4.6. Tetragonal symmetry 

resentation: 

/ 



D 



\ 

where {5, a) G M?. The monodromy group is defined by: 

6^5, (7 I— )• —a. 
The evaluation of the invariants (16) on this slice gives: 
J2 = 8 280 (5^ 



a + 35 


a + 5 


-4 6 











a + 6 


-a + 36 


-4 6 











-4 6 


-45 


86 




















-86 




















-8 6 




















2a + 26 



\ 



(25) 



J3 

J5 
Je 
J7 
Js 
J9 
Jio 

We notice that the parameter 



48 6 {a^ + 15 6^) 
32 0-^-^960(5^2-^ 32800(5' 
128(5 {56-af{56 + af 



128 0-^ + 5760 0-^(5^ + 86400 a^6'^ + 4528000 6^ 
512 (5 (55 6^ + cj2) {5 6- af (5 5 + cj)^ 
6144 (5^ (55 (5^ + cj2) (5 (5 - (7)^ (5 (5 + (j)^ 
2048 6 {b6-af{b6 + of (55 (5^ + a^f 
24576 (52 (5 ,5 - (j)^ (5 (5 + cj)^ (55 (5^ + a'^f 



1 



4 J22 - 3 J4 
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is a rational invariant, therefore for all {J2, Ji, J5) G 
satisfied by a is of degree 2 
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However, the minimal equation 



J2 = 8 + 280 5^ 



As for [B3], the condition 



2J2(J| -3J4) 



2 35J| > 



(26) 
(27) 



has to be added in order to ensure a to be real. 

Having computed 6, the equation (26) has two roots with opposite signs. The linear slice corre- 
sponding to the group D4 is of degree 2, which was already known (see Appendix B). 



Proposition 4.9. A harmonic tensor D £ M.^ is in the closed stratum S[ 
invariants J2{D) ■ ■ ■ Jio{D) satisfy the following polynomial relations 



if and only if the 



6J6 

3 J7 

6 

5J8 

3 

6J9 
5J5' 
15Jio 



-3J2'' + 9J2J4 + 20 Js^ 

J2^Jz — 3 J^Ja + 3 J2 J5 
-3 J2^ + 9 J2 V4 + 20 J2 Ja^ - 20 J3 J5 
-3 J2^ + 6 J2 V4 + 30 J2 Ja^ - 5 J3 J5 
15 J2J3J4 + 3 J2V5 
I3J2J3J4 + 6 J2V5 
h 



7J2V3 
5J2V3 



60 Js^ 
20 J,^ 



(28) 



2 J2^ + 4 J2 V4 + 10 J2 - 20 J2 J3 J5 + 30 J3 V4 



-9 J2^ + 18 J2 V4 + 85 J2 V32 - 30 J2 J3 J5 + 15 J3 V4 



together with inequality (27) . D is in the tetragonal class [B4] z/ moreover 3 J4 — Jl 7^ antf 
5 JI — 8 J2J4 — 70 J| / 0. /n i/iai case, it admits the normal form (25) where 



1 



^5 



4 Js^ - 3 J4 

and where a is the positive root of J2 = 8 fi^ + 280 5^ . 

Remark 4.10. All 2nd-order tensors in covariant to D are multiple of diag(l,l,— 2) and are 
eigenvectors of D corresponding to the eigenvalue 125. 



4.7. Orthotropic symmetry ([ID'2])' A 
representation: 



D 



-invariant tensor D £ has the following matrix 





— A2 — A3 


A3 


A2 








\ 




A3 


—A3 — Ai 


Ai 













A2 


Ai 


— Ai — A2 






















2Ai 






















2A2 





V 

















2A3 / 



(29) 



where (Ai, A2, A3) G M'^. The monodromy group is the group 63 acting on (Ai, A2, A3) by permuta- 
tion. The evaluations of the invariants (16) on this slice are ©s-invariant. They can therefore be 
expressed as polynomial functions of the elementary symmetric polynomials (o"i, (T2, (73) defined by: 



fji := Ai -|- A2 + A3 



(T2 :— A1A2 + A1A3 -I- A2A3, (T3 :— A1A2A3. 
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More precisely, we have: 

J2 = -14(72 + 8(7? 
J3 = —6(71(72 + 24(73 

J4 = 40 (71(73 - 112 (7^(72 + 68 cri + 32 aj 

J5 = 64 (7i (73 - 12 (72(73 - 16 (7^(72 + 28 (71(7^ 

Je = -344 (7^ + 192 (7? (73 - 24 erf - 672(7^(72 + 1008 erf (7^ + 128 af - 504 (71(72(73 

J7 = -432 (7^(72(73 + 384 afas + 104 (7^(73 - 96 aiaj - 64 (7^(72 + 192 afa^ - 248 aia^ 

Js = 608 (7? (72(73 + 80(7^ - 768(7^(73 + 192 afa'^ + 72(72(7| + 288 erf erf - 416(7^(71 + 744(71(7^(73 

Jg = -5248(7^(72(73 + 28 80 (7? (7^ (73 + 1328 (7l(72(7| + 144 (7| + 2304(7^(73 - 1152 (7? erf - 880(7^(73 

- 256(7^(72 + 1024 (7^ (7^ - 2304 (7? (7^ + 2160 (7i (7^ 
JlO = 10752 (7i (72(73 — 1280(7^(7^(73 — 5664(71(7^(73 — 2688 (7^ (72 crf — 800 (7| -4608(7^(73 
+ 2304(7f(7| - 1344(7|(7| - 288(7i(7| + 1536 (7? (7^ - A224afa^ + 3104 (7? (7^ 

As it can be observed, the parametric system obtained for [D2] is rather comphcated. Nevertheless, 
the evaluation of J2,..., Jio on the slice and their expressions using ai, (72, (73 is simple using 
computer. By computing the elimination ideal of a Groebner basis for the ideal generated by 

Jk-Pk{<^i, 0-2,0-3), k = 2,...,7, 

we obtain a set of 6 syzygies. 

- 1350 J3 J7 - 840 J4 Je + 465 J2^ Je + 270 J^"^ + 720 J2 J3 J5 + 747 J2 

- 170 J32 J4 - 564 J4 + 70 J2^ J3^ + 84 J2^ = (30a) 

- 1620 J4 J7 + 810 J2^ J7 + 360 J5 Je - 1110 J2 J3 Je + 999 J2 J4 J5 + 960 Ja^ J5 

- 549 J2^ J5 - 972 J3 J42 + 1638 J2^ J3 J4 - 80 J2 J3^ - 312 Ja'' J3 = (30b) 

4050 J5 J7 - 25650 J2 J3 J7 - 14310 J2 J4 Je + 9600 Ja^ Je + 7965 J2^ Je + 9450 J3 J4 J5 
+ 10530 J2^ J3 J5 + 1134 J4^ + 11259 J2^ Ja^ - 12330 J2 J3^ J4 - 9018 Ja^ J4 

+ 400 Js"^ + 3270 J2^ J3^ + 1350 J2^ = (30c) 

- 12150 J2 J3 J7 + 3600 Je^ - 11610 J2 J4 Je + 9750 J3^ Je + 4410 Ja^ Je 

+ 8505 J3 J4 J5 + 3645 J2^ J3 J5 + 1458 J4^ + 5670 Ja^ Ja^ 

- 10710 J2 J3^ J4 - 4104 J2^ J4 + 400 Ja"^ + 2580 Ja^ J3^ + 576 J2^ = (30d) 



1800 Je J7 - 10800 J2 J4 J7 + 4800 J3^ J7 + 4950 J2^ J7 + 4020 J3 J4 Je - 8370 Ja^ J3 Je 
+ 162 J4^ J5 + 7371 J2^ J4 J5 + 2880 J2 J3^ J5 - 3483 J2'' J5 - 9216 J2 J3 J4^ 

+ 640 J3^ J4 + 11946 J2^ J3 J4 - 720 J2^ J3^ - 2160 Ja^ J3 = (30e) 
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60750 Jr^ + 178200 J3 J4 J7 - 546750 J3 J7 + 3780 Ja^ J& - 246780 J4 Je 

+ 348000 J2 Js^ Je + 137025 Je + 116640 J2 J3 J4 J5 - 75600 Jg^ J5 

+ 223560 J2^ J3 J5 + 29808 J2 J4^ + 82170 Ja^ J4^ + 177660 J2^ Ja^ 

- 438390 J2^ Ja^ J4 - 148014 J2^ J4 + 17200 J2 Ja^ + 102000 J2^ Ja^ + 22221 J2'' = (30f) 

If these relations are satisfied, the elementary symmetric polynomials ui, (T2 and ua can be 
expressed, on the open stratum defined by 

6 Je - 9 J2 J4 - 20 Ja^ + 3 J2 V 0, 
as rational expressions of J2, . . . , J7. More precisely, we have 

4 2 1 r 1 r 1 3 1 

a2 = -ai --J2, a3 = -Ja + -ai -^^1^2, 

and 

_ 9 (3 J7 - 3 J2 J5 + 3 Ja J4 - J2^ Ja) 
2 (6 Je - 9 J2 J4 - 20 Ja^ + 3 J2^) ' 

Then Ai,A2,Aa are recovered (up to their monodromy group) as the roots of the third degree 
polynomial in A 

p{X) = A^ - ai\^ + cr2A - da 

But to ensure that given J2, . . . , J7 G M, we can find Ai, A2, Aa € M, the following inequalities have 
to be satisfied 

Aa = alal - 4cj?cra + 18cJicJ2cr| - 27cr| > , 

as well as 

A2 = 2al - 60-2 > . 
The discriminant Aa can be rewritten, using J2, ■ ■ ■ ,Je, as 

6 Je - 9 J2 J4 - 20 Ja^ + 3 Ja^ 

ZAa = . 

432 

whereas A2 can be rewritten, using J2, . . . , J7, as 

iV2 



A. 



14 (6 Je - 9 J2 J4 - 20 Ja^ + 3 J2^f 



where 



N2 = 3645 J72 + 7290 J2 J5 J7 - 7290 Ja J4 J7 + 2430 J2^ Ja J7 - 36 J2 Je^ 

+ 108 J2^ J4 Je + 240 J2 Ja^ Je - 36 J2'' Je - 3645 J2^ Js^ + 7290 J2 Ja J4 J5 
- 2430 J2^ Ja J5 - 3645 Ja^ J4^ - 81 J2^ J4^ + 2070 J2^ Ja^ J4 
+ 54 J2^ J4 - 400 J2 Ja^ - 285 J2^ Ja^ - 9 J2^. 
We will summarize these results in the following proposition. 

Proposition 4.11. A harmonic tensor D E is in the closed stratum Sp^] if and only if the 
invariants J2{D), ■ ■ ■ , Jiq{D) satisfy the polynomial relations (30), completed with the 3 equations 

5fc(J2,... Jio) = 0, A; = 8, 9, 10 (31) 

obtained from J^ — Pfe(o"i, fT2, ua), where we have substituted the fij by their corresponding rational 
expression, and we have cleared the denominators. The inequalities 

6 Je - 9 J2 J4 - 20 Ja^ + 3 J2^ > (32) 
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as well as 

3645 Jt^ + 7290 J2 J5 J7 - 7290 J3 J4 J7 + 2430 J3 J7 - 36 J2 Je^ + 108 J4 Je 
+ 240 J2 Js^ Je - 36 72"^ Je - 3645 J2^ Js^ + 7290 J2 J3 J4 J5 - 2430 J2^ J3 J5 

- 3645 Ja^ J4^ - 81 J2^ J4^ + 2070 J2^ Js^ J4 + 54 Js^ J4 ^^^^ 

- 400 J2 Js"^ - 285 72"^ Ja^ - 9 J2^ > 

are ako required. D is in the orthotropic class [ID2] if moreover inequality (32) is strict. In that 
case, it admits the normal form (29) where Ai,A2,A3 are the roots of the third degree polynomial 

p{X) = - aiX^ + a2X - 0-3. 

Remark 4.12. Even if it is the core of the method, it worths to note that, instead of Ai, A2, A3, the 
3 numbers o"i , fT2 , CT3 are uniquely defined on the open stratum: they are invariants of D and can be 
written as rational expressions of {J2, ■ ■ ■ , Jy). It should be noticed that it is also possible to express 
ci, <T2, <T3 as rational functions of only {J2, ■ ■ ■ , J5) on an open subset of Sp^j. We will, however, not 
give these formulas here. 

Remark 4.13. All second rank tensors in covariant to D commute together and the commutators 
of two of any of them vanishes. 

4.8. Bifurcation conditions for tensors in H^. An important observation has to be made. 
Mathematical results [34] tell us that the invariant algebra of is generated by the 9 fundamental 
invariants ( J2, . . . , Jio)- In this set 6 invariants ( J2, . . . , J7) are algebraically independent, meanwhile 
the others are linked to the formers by polynomial relations. As discussed before, it is a well-known 
fact that this algebra separates the orbits. But as demonstrated exhaustively in this section, for any 
class with a finite monodromy group^*', the first 6 invariants (J2, . . . , J7) are necessary and sufficient 
to separate the orbits inside each class. This observation can be summed-up by the following 
important theorem. 

Theorem 4.14. For the following isotropy classes: the cubic, the transversely isotropic, the trig- 
onal the tetragonal (and the trivial isotropic^ classes in IH[^(R^), the first 4 invariants (J2, . . . , J5) 
separate the orbits inside the class. For the orthotropic class, the first 6 invariants (J2,..., J7) 
separate the orbits inside the class. 

We will now summarized the previous results by giving bifurcation conditions on the to explicit 
how to "travel" from a given isotropy class to another. More precisely, so far, we have given necessary 
and sufficient conditions to belong to a closed stratum "^[h] with finite monodromy group. As we 
have seen, ^[h] is characterized by a set of polynomial relations (syzygies) involving the invariants 
Jfc. These relations define implicit equations which characterize the closed stratum The set of 
relations that characterizes a given (closed) stratum generates an ideal (c.f. subsection A. 4), I^ijt^ 
in M[J2, . . . , Jio], and we have 

The bifurcation conditions, detailed on figure 2, correspond to relations (generators) which need to 
be added to I^jj-^j to obtain I[H2] ■ 

Remark 4.15. These bifurcation relations were obtained with the following algorithm. Starting with 
some subgroup H and the corresponding fixed point space , we look for linear relations needed 



i.e. for all symmetry classes, except the monoclinic and the triclinic classes 
^-'^Notice that a rigorous and complete classification of these closed strata would require also some real inequalities 
that we have not completely detailed here. 
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(T) isotropic: [S0(3)] 




Figure 2. Bifurcation paths for finite monodromy isotropy classes in M.^. 



to define the subspace of , where H C K. Given such a relation a = (in fact, in all the 
cases considered, only one relation was needed), we compute a F-'^-invariant polynomial 

provided F''^ is finite. It is then clear that pa vanishes on a point v G if and only if v belongs to 
or yT^T , for some 7 G F^, that is if and only if v has at least isotropy [K]. This F^-invariant 
polynomial is then expressed using the Jfc. 

4.9. Relation with the characteristic polynomial. To conclude this subsection, some links 
that exists between the invariants of H'* and the characteristic polynomial of D will be investigated. 
The first link is formulated by the following proposition: 

Proposition 4.16. Let D, G H^. The coefficients of the characteristic polynomial of D, are S0(3)- 
invariants of D. They can be expressed in terms of the invariants J2,...,J5 defined above. We 
have 

Xd{z) = - ^ J2Z^ - I J3Z^ + I {Ji" - 2 J4) z2 + A ( J2 J3 _ 3 J5) z. (34) 

Proof. Notice first that 2: = is always an eigenvalue of D^. Indeed, if q stands for the standard 
metric on M^, we have Z2q = tri2(-D) = 0. The coefficient of in XD is a homogeneous invariant 
polynomial of degree 6 — A; for 1 < /c < 4. Therefore, it can be expressed as a polynomial in 
J2{D), ■ ■ ■ J^{D). The computation of this last expression reduces therefore to identify some real 
coefficients. This can be done by specializing on particular values or on some particular strata where 
it is easy to compute, the cubic or the transverse isotropic stratum, for instance. □ 

Besides the characteristic polynomial, Betten [8] has introduced a two variables, invariant poly- 
nomial of the elasticity tensor. It is defined as 

Sc(A,/i) = det(C-C,,^) 
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where 





A + 2/i 
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A + 2/X 
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A + 2;U 
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2fi 





v 

















2^ / 



is the totahy isotropic tensor. Notice that ;Bd(0, /u) = xd(2//) and that Be is of degree lower than 
1 in A (indeed, this can be observed by subtracting the first hne to the second hne and to the third 
hne in the determinant). The supposed interest of this polynomial is to obtain more invariants 
than the ones of the characteristic polynomial. However, it seems that it does not bring any new 
invariants for a tensor D £ IH^(M^). Indeed, we have the following result. 

Lemma 4.17. The Betten polynomial of a harmonic tensor D G E[^(M'^) is related to its charac- 
teristic polynomial by the formula 

BD{X,f^) = (3A + 2/X) xi)(2^) 

where xd(A) = XxdW- 

Proof. Let {I < k < 6) denote the lines of the determinant det{D_ — C_x^^)- For a harmonic tensor 
D, we have 

3 

^ dp,fc = 0, I <p< 6. 

k=l 

Thus, substituting Li — )• Li + L2 + ^3 in the determinant det(^ — C_\^^), we get 

BD{X,^^) = (3A + 2/i)P(A,/i) 

but since Bd{X, /x) is of degree lower than 1 in A we have -P(A, //) = P(0, ^). Now for A = we have 

^d(0,/x) = Xd(2/x) = (2^)P(0,;u) 

from which we get P{0,^) = x'^{2fj,). This achieves the proof. □ 

Corollary 4.18. The invariants defined by the coefficients of the Betten polynomial do not separate 
the orbits o/Ela. 

Proof. If it was the case, it would also separate the orbits of but then lemma 4.17 and proposi- 
tion 4.16 would imply that J2, J3, J4, J5 separate the orbits of H^, which is false. □ 

5. Towards an effective determinacy of elasticity classes 

In this last section, an application of the concepts previously introduced will be given. This 
application will concern the identification of the symmetry class of a tensor which is not expressed 
in one of its natural bases. As our objective is mainly to illustrate the tools we propose, we will 
only provide a set of sufficient conditions. This approach is based on a hypothesis of genericity for 
the anisotropic tensors. Such hypothesis, which amounts to exclude any exceptional case, will be 
detailed and criticized in the next subsection. Therefore, in order to avoid any misunderstanding, 
our aim is not to provide a complete answer to identification problem. A full treatment of this 
question would require a detailed analysis of each particular non generic case that may occur, and 
constitutes a full and complex study on its own. As the core of this first paper is focused on the 
symmetry class identification in terms of EI'* invariants, we have rather choose a simpler problem to 
treat. 
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5.1. Generic anisotropic elasticity tensors. A well-known result of anisotropic elasticity [19] 
states that the symmetry group of an elasticity tensor is the intersection of the symmetry group of 
each one of its irreducible components. In other terms: 

Gc = Ga^Ghn Gd 

since we trivially have Gx = Gfj_ = S0(3). In the most general situation, all the former tensors 
are independent and therefore a complete analysis of the intersection of symmetry group has to 
be conducted. Such a problem is a complex one, and is a subject of a study on its own. For the 
physical situations considered here, a simplifying hypothesis will be made. As the tensor symmetry 
is a consequence of the material symmetries, a weak coupling between the tensors of the set (a, b, D) 
will be considered. To that end let us consider the triplet of quadratic forms L := (a, 6, d2), where 
d2 is the fundamental covariant of degree 2 introduced in (16). We have the following direct lemma 

Lemma 5.1. 

Gc = Gl n Gd 

Proof. As d2 is covariant to D, G^) C and therefore n Gjj = Gd- Thus 
Gl n G/) = Ga n G& n Gd2 n G/) = Ga n G^ n Gd = Gc 

□ 

The hypothesis that will be made here is the following 

Hypothesis 5.2 (Material Genericity Assumption (MGA)). // [Gl] >z [D2] then Gl = Gda 

Such an hypothesis is equivalent to the following: 

if [Gl] h [D2] then Gd^ C Ga and Gd^ Q G^. 

From a physical point of view, this means that the symmetries of a and h are compatible with the 
ones of D in the sense that adding Ga and G^, to the intersection will not lower Gd- In other words 

Lemma 5.3. 

Gl = Gda =^ Gc = Gd 

Proof. From the lemma 5.1, we have Gc = Gl H Gd- Therefore, according to the hypothesis 
Gc = Gd2 n Gz3, and as d2 is covariant to D the intersection reduces to Gc = Gd- D 

If this genericity assumption fails, and without a deeper study, all it can be said is Gc ^ Gd- 
From a mathematical point of view, studying these cases are necessary to close the problem. On the 
other hand, from a material point of view and for a first approach, we believe that the genericity 
assumption to hold true for a broad class of classical materials. Anyway situations where this 
hypothesis may fail can be conceived: 

Tailored materials: As the symmetry class of an elasticity tensor is driven by the architec- 
ture of the microstructure, it can be conceived to tailor the microstructure in order to give 
the material some non-conventional properties [25]. Therefore for such a specific design the 
genericity can not be assumed. Such kind of very specific architecture might not only be 
found in engineered materials, but also in some biomaterials. 

Damaged materials: Genericity may also fail for damaged-materials, according to certain 
model [22] the symmetry breaking induces by the damage process may result in non generic 
anisotropic properties. 

Therefore, bearing in mind the aforementioned limitations induced by the genericity assumption, 
we can propose, in the following section, an identification algorithm. 
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[Gl] = [1] ? 



[Go] = [1] 



C 

\ 

{a,h,D) 
i 

[Gl] = [Z2] ? [Gl] = [B2] ? 



[Gc] = [1 



[Gc] = [1] 



[Gl\ = [0(2)] ? [Gl] = [S0(3)] ? 

- [Go] = [0(2)] ? U [Gd] = [S0(3)] ? 
i i 
[Gc] = [0(2)] I [Gc] = [80(3)1 



Yes 



[C?d] = [O2] '! 



.,1 ? No 



[Gc] = [IO2] 



[Gc]e{[l].[Z2]} 



[Go] = [D4] ^ 



■ [Go] = [O]? 



[Gc] = [D4: 
Go] = [B3] 



Gc] = [O] 
[Go] = [02] ? 



[Gc] = 



[Gc] = 



[Go] = [B2] 

I 

[Gc] = [»2] 



else 



[Gc] S{[1],[Z2]} 



else 
i 



[Gc]e{[t],[U] 



Figure 3. Symmetry identification for a generic tensor in Ela. 



5.2. Symmetry class identification. Based on the results obtained in tliis paper, tfie algorithm 
provided figure 3 allows to identify the symmetry class under the MGA assumption. Concerning 
this algorithm some remarks are in order. 

• The first step of the algorithm is based on the symmetry analysis of the triplet of quadratic 
forms L. This symmetry identification is based on the conditions given in subsection 2.2. As 
indicated on the diagram if Gl = ^2, the symmetry class of C is decided using the Cowin- 
Mehrabadi criterion for the non vanishing vector cj, obtained from case (1) of Theorem 2.1 
applied to the triplet L. This test is denoted by CM(a;)? in figure 3. 

• For classes strictly greater than [D2] we use the syzygy relations given in section 4, to check 
the isotropy class of D. 

Some tests can be avoided a priori (and will thus not be included in figure 3) by the following 
observations which result from lemma 5.4 combined with the bifurcation conditions given in figure 2. 

(1) If Gl = 0(2) = Gd2, then the isotropy class of D cannot belong to {[O], [S0(3)]}. 

(2) If Gl = S0(3) = Gd2: then j| — 3J4 = and the isotropy class of D cannot belong to 
{[D3],[B4],[0(2)]}. 



Lemma 5.4. //d2 



Aid then j| 



3J4 = 0. 



Proof. If d2 is a scalar multiple of the identity, Jl 
conclusion follows. 



9A2 and J4 = tr(d4) = tr(d^) = ?>\^ 



The 

□ 



6. Conclusion 



In this paper a method aiming at identifying the symmetry class of an elasticity tensor has been 
provided. Conversely to other approaches based on the spectral decomposition of the elasticity 
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tensor, this method is based on the specific relations between invariants in each symmetry class. A 
complete set of relations were given for the couple of quadratic forms appearing in the harmonic 
decomposition of Ela, and 5 sets of algebraic relations were provided to identify the following 
symmetry classes from H^: orthotropic ([ID)2])) trigonal ([lUs]), tetragonal ([1154]), transverse isotropic 
([S0(2)]) and cubic ([O]). Completed by a genericity assumption on the elasticity tensor under 
study this approach was shown to be able to identify the 8 possible elasticity classes. 

We consider this paper as a first step towards a more systematic uses of invariant-based methods 
in continuum mechanics. Using the geometrical framework introduced in this paper some extensions 
of the current method can be envisaged : 

• The first point would be to get rid of the genericity assumption we made, such a point 
constitutes a full study on its own ; 

• We aim at extending the current sygyzy-characterization to group having continuous mon- 
odromy group ; 

• To propose a complete set of normal forms of elasticity tensor constructed using this invari- 
ants approach ; 

• To tackle the increasing degree of the algebraic relations for low symmetry group, it might 
be interesting to characterize classes in terms of covariants. 

Such an approach was conducted here in the context of classical elasticity, but the framework 
used allows a direct extension of this method to other kinds of anisotropic tensorial behaviors. 
It can be suited for the study of anisotropic features of piezoelasticity [23], flexoelasticity [30] or 
strain-gradient elasticity [29]. Furthermore, and in a more experimental view, we aim at testing our 
approach on experimental situations. This would allow to compare the proposed approach with the 
more classical ones. 

Appendix A. Geometric concepts 

A.l. Normal form. A normal form is a map NF : V ^ V such that 

vi w V2 if and only if NF(wi) = NF(t;2). 

The existence of a computable normal form allows to decide whether two vectors vi and V2 belong 
to the same G-orbit. 

Remark A.l. In practice, one contents itself usually with a normal form defined up to a finite 
group. For instance, one is satisfied with a diagonal matrix as "the" normal form for the action of 
the rotation group on the space of symmetric matrices. However this diagonal form is not unique. 
It is defined up to the permutation of the diagonal entries. 

A. 2. Covariants. Given a second representation {V', p') of the same group G, we say that a linear 
map if '.V ^ V' is equivariant if the following diagram 



V v 



pia) 



p'ia) 



V V 

commutes for each g £ G. The composition of two equivariant maps is equivariant. 

Given v € V the element ip{v) will be said to be covariant to v. It is clear that every element 
V £ V is covariant to itself (with ip = Id). However, this relation is not symmetric. Indeed, if v' is 
covariant to v, the converse is not necessarily true {v is not necessarily covariant to v') unless the 
linear map if : V — > V' is bijective. In this case, we say that the representations {V, p) and iV' , p') 
are isomorphic or equivalent. 
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When the representation {V',p') is trivial, i.e. when p'{g) = Id for ah g £ G, the covariant 
element '^{v) wih be cahed an invariant. We have then 

V M v' ^ ip{v) = (p{v'). 

A. 3. Polynomial invariants. The hnear action of the group G on the M-vector space V extends 
naturahy to the vector space of polynomial functions defined on V with values in M. This extension 
is given by 

ig-P){v) ■.= P{g-^.v) 

for every polynomial P £ M[V] and every vector v £ V. The set of all invariant polynomials is a 
sub-algebra of the algebra M[V] of all polynomials defined on V. We denote this sub- algebra by 
M[y]'^ and call it the invariant algebra. 

When G is a compact Lie group (or more generally a reductive group), this algebra M[l/]*-' is 
of finite type [27]. This means first that it is possible to find a finite set of invariant polynomials 
Ji, . . . , Jtv which generate (as an algebra) R[l/]'^. But, in general, this algebra is not free, which 
means that it is not always possible to extract from Ji , . . . , a subset of generators which are 
algebraically independent^^. The set of polynomial relations on Ji, . . . , Jn (also known as syzygies) 
is an ideal of M[Ji, . . . , J^r] that is finitely generated. Such a set of generators is called an integrity 
basis for ]R[y]'^. Of course it is not unique and moreover its cardinal may change from one basis to 
another. 

A. 4. Ideal and Groebner basis. Consider a subset V of A;" (where = M or C) defined implicitly 
by a finite set of polynomial equations pi, . . .pr- The fundamental observation is that the set / of 
polynomials p £ k[xi, . . . ,Xn] such that p{x) = for each point x £ S (where x = (xi, . . . ,Xn)) 
satisfies the following properties: 

(i) £ I, 

(ii) li p,q £ I, then p + q £ I, 

(iii) li p £ I and q £ k[xi, . . . , then pq £ I. 

Such a subset / of k[xi, . . . , Xn] (or more generally any commutative ring) is called an ideal. In the 
example above the ideal / is (finitely) generated hy pi, ... pr (i.e. every p £ I can be written as a 
finite sum X]j=i fjPj where fj £ k[xi, . . . , Xn]) and we write I =< pi, . . .pr >■ 

To solve a system of linear equations in k", a standard algorithm is to reduce the matrix of the 
system to a triangular form (a normal form, even if not unique). A Groebner basis of an ideal / is 
a set of generators for / which is a normal form for an algebraic system, in a certain sense analog 
to a triangular system for a set of linear equations. It permits to compute and answer eff'ectively 
a lot of problems related to ideals k[xi, . . . ,Xn] and algebraic systems. It solves, in particular, the 
implicitization problem. We refer to [17] for an excellent exposition on Groebner basis (which is 
very accessible for non specialists). 

A. 5. Separation of orbits. Given an algebra A of invariant functions defined on V, we say that 
A separates the orbits if given two vectors vi,V2 £ which are not on the same orbit, it is 
always possible to find a function f £ A such that f{vi) 7^ f{v2)- In general, when the group G 
is not compact, the algebra of polynomial invariants does not separate the orbits^^. However, it is 
remarkable that for any compact group, the algebra of real invariant polynomials always separates 
the orbits. This result is attributed to Weyl [37]. 



Polynomials Pi , . . . , Pjv are said to be algebraically independent over R if the only polynomial in A'^ variables 
which satisfies Q{Pi, . . . , Pjv) = is the zero polynomial. 

-'^'^This can be illustrated by the example of the adjoint action of the general linear group GL(n, C) on its Lie 
algebra 0l(n, C). Each orbit is described by a unique Jordan form but the algebra of polynomial invariants, which is 
generated by tr(M), tr(M^), . . . tr(M") where M £ £|l(n, C), is not able to distinguish between two different Jordan 
forms which have the same diagonal part. 
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Proposition A. 2. Let (V, p) be a finite dimensional, linear representation of a compact group G. 
Then the algebra of real invariant polynomials separates the orbits. 

Proof Let v €z V and K a compact subset of V. Given a norm ||-|| on V, we define the distance of 
f to as 

d{v, K) = min||f — 

Now, let Ki and K2 be two distinct G-orbits in V . The following function 

^ d{v,Ki)-d{v,K2) 
d{v,Ki) + d{v,K2) 

is continuous on V , f{v) = —1, for v £ Ki and f{v) = +1, for v £ K2. 

Since the group G is compact, the orbits are compact and it is possible to find a closed ball B 
in V which contains both orbits. We apply then the Stone-Weierstrass theorem which states that 
real polynomial functions are dense in the space of continuous functions defined on B (for the sup 
norm). Therefore, a polynomial function PonV can be found such that 

P{v) < e Ki, P{v) > ^,yv e K2. 

Now, averaging this polynomial P 

P= [ g-Pdii 
Jg 

where is the Haar measure over G, we obtain a G-invariant polynomial P which separates the 
orbits Ki and i^2- D D 

Appendix B. Normalizers of closed subgroups of S0(3) 
Every closed subgroup of S0(3) is conjugate to one of the following list [18, 24] 

S0(3), 0(2), S0(2), D„(n > 2), Z„(n > 2), T, O, I, andl (35) 

where: 

• 0(2) is the subgroup generated by all the rotations around the z-axis and the order 2 rotation 
a : (x, y, z) 1— )• (x, — y, —z) around the x-axis. 

• S0(2) is the subgroup of all the rotations around the z-axis. 

• Z„ is the unique cyclic subgroup of order n of S0(2), the subgroup of rotations around the 
z-axis. 

• D„ is the dihedral group. It is generated by Z„ and a : (x, y, z) 1— )■ (x, — y, —z). 

• T is the tetrahedral group, the (orientation-preserving) symmetry group of a tetrahedron. 
It has order 12. 

• O is the octahedral group, the (orientation-preserving) symmetry group of a cube or octa- 
hedron. It has order 24. 

• I is the icosahedral group, the (orientation-preserving) symmetry group of a icosahedra or 
dodecahedron. It has order 60. 

• 1 is the trivial subgroup, containing only the unit element. 

The lattice of conjugacy classes of closed subgroups of S0(3) is completely described [24] by the 
following inclusion of subgroups (which generates order relations between conjugacy classes) 

Z„CD„C0(2) (n>2), 

Z„ C Zm and ID)„ C B^, (if n divides m), 

Z2 CBn (n > 2), 

Zn C S0(2) C 0(2) (n > 2), 
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and the arrows in figure 4 which complete the lattice, taking account of the exceptional subgroups 
0,T,I (beware that figure 4 cannot be realized by an inclusion diagram between representatives of 
these conjugacy classes). 



[S0(3)] 




[D4] [T] [D3] [B5] 




[1] 



Figure 4. Exceptional classes in the lattice of conjugacy classes of closed subgroup 
of S0(3). 



Remark B.l. Notice that O and I are maximal subgroups of S0(3) which both contained T. 

Recall that the normalizer of a subgroup H of a group G, denoted by N{H) is the biggest 
subgroup of G in which H is normal. It is defined by 

N{H) := {g G G; gHg-' = H] . 

Proposition B.2. The normalizers of the subgroups in the list (35) are given below. 

Ar(S0(3)) = S0(3), N{0{2)) = 0{2), A^(S0(2)) = 0(2), 

N{0) = O, N{I) = I, N{T) = O, 

A^(1) = S0(3), N{Zn) = 0{2) for n> 2, iV(B„) = Ban /orn > 3, 

and N{I}2) = O. Moreover, the quotient groups are given by 

7V(SO(3))/SO(3) = 1, iV(0(2))/0(2) = 1, iV(SO(2))/SO(2) = Z2, 

iV(0)/0 = l, N{I)/I = 1, N{T)/T = Z2, 

N{1)/1 = S0(3), iV(Z„)/Z„ = 0(2) for n > 2, iV(B„)/D„ = Z2 for n > 3, 

and N{I}2)/02 = ©3, the symmetric group of 3 elements. 

Proof. That iV(S0(3)) = S0(3) and N{\) = S0(3) are obvious. The normalizers of 0(2), O and I 
are equal to themselves because each of these groups is a maximal subgroup of S0(3) and because 
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S0(3) is a simple group (it has no non-trivial normal subgroups). Given any rotation G S0(2), 
we have 

area = r_g 

where a : (x,y,z) i— )■ (x, —y, —z). Hence 0(2) c A^(S0(2)) but since 0(2) is maximal and S0(3) is 
simple we get that A''(S0(2)) = 0(2). T is a subgroup of index 2 of O and is thus a normal subgroup 
of O. Therefore, we have O C N{T) but since O is maximal and S0(3) is simple, this proves that 
A''(T) = O. Next consider Z„ for some n > 2. Since Z„ is a subgroup of the abelian group S0(2) 
we have S0(2) C N{Zn). Moreover 

0"''27r/nO" = r_27r/n ^ Z„ 

and therefore 0(2) C N{Zn). Using again the maximality argument, we get N{Zn) = 0(2). Finally, 
consider D„ for some n > 2. Because D„ is subgroup of index 2 of 02n, we have D2ra C A^(D„). 
Suppose first that n > 2 and let g £ A^(D„). Then gr2.,r/n9~^ is a rotation of angle =b27r/n and its 
axis is therefore the 2;-axis. Hence g{es) = ±63. Besides, gcrg~^ is a rotation of order 2 whose axis 
is not the z-axis (otherwise we would have g{ei) = ±63). Therefore we have 

for some k G {0, 1, . . . , n — 1}. Now if g{e^) = 63, we have g = r^ and necessarily 29 = Akir/n. If 
5(63) = —63 we have g = are and necessarily 26 = Akir/n. In both cases, this means that g G 
and we conclude that Ar(D„) = D2n- Suppose now that n = 2, in which case D2 is the Klein group 

^2 {e, Rei , Re2 ! Re.i } , 

where Rg. denotes the rotation by angle tt around ej. For each element g of the octahedral group 
O, we have g{ei) = ztcj and hence 

9Rei9~^ = Rey 

Therefore O C A''(B2). Using again the maximality argument, we conclude that N{J])2) = O- 

For the computation of the quotient groups, all of them are obvious but Ar(Z„)/Z„ and N{p>2)/'^2- 
Concerning the first quotient group, we remark that every element in 0(2) can be written uniquely 
as 

where £ = 0, 1 and r G S0(2) is a rotation around the z-axis. Consider now the map 

(j) : 0(2) ^ 0(2), a^r ^ a^r"". 

One can check that this map is a surjective, group morphism and that its kernel is exactly Z„. This 
establish the isomorphism between 0(2)/Z„ and 0(2). 

For the second quotient group, we recall that there is a well-known isomorphism from the sym- 
metry group of the cube O onto the group of permutations of the diagonals of the cube, that is 64. 
Under this isomorphism, the subgroup 1^2 of O is sent to 

V2 := {e, (12) (34), (13) (24), (23) (14)} , 

which is a normal subgroup of ©4. Let ©3 be the subgroup of ©4 which fixes the element 4. Then 
©4 is the semi-direct product of T>2 by ©3. In particular, the quotient ©4/X'2 is isomorphic to ©3. 
This subgroup ©3 corresponds geometrically to the subgroup of O generated by a rotation of angle 
27r/3 around a diagonal and a rotation of angle tt around any diagonal perpendicular to the first 
one, that is a conjugate of the group D3. □ 
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Appendix C. Trace formulas for 



Let (y, p) be a finite dimensional linear representation of S0(3), which admits the following 
harmonic isotypic decomposition : 



(36) 



fc=0 



where is the multiplicity of the space M.^ in the decomposition. Using the trace formula (3) the 
following explicit formulas can easily be obtained [24, 4] 
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